plot 1-(R^2/2 + R + 1) *e^(-R) with R from 0 to 10